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Abstract. We study C 2 extremal quasiconformal mappings in space and establish neces- 
sary and sufficient conditions for a 'localized' form of extremality in the spirit of the work 
| of G. Aronsson on absolutely minimizing Lipschitz extensions. We also prove short time 

existence for smooth solutions of a gradient flow of QC diffeomorphisms associated to the 
| /" ( | extremal problem. 

JT" 1, 1. Introduction 

Ph ; 

A quasiconformal (qc) mapping is a homeomorphism, u : f2 C W 1 — > W 1 whose components 
are in the Sobolev space and such that there exists a constant K > ^fn for which 

\du\ n < Kdetdu a.e. in Q. Here we denote |v4| 2 = Y17j=i a ij ^° ^ e the Hilbert-Schmidt 
norm of a matrix and du the Jacobian matrix of u = (it 1 , ...,u n ) with entries duij = djU 1 . 
At a point of differentiability du(x) maps spheres into ellipsoids and the smallest possible K 
in the inequality above, roughly provides a bound for the ratio of the largest and smallest 
axes of such ellipsoids. In this sense qc mappings distort the geometry of the ambient space 
in a controlled fashion. Quoting F. Gehring [19] . qc mappings " constitute a closed class 
of mappings interpolating between homeomorphisms and diffeomorphisms for which many 
results of geometric topology hold regardless of dimension." 

Quasiconformality can be measured in terms of several dilation functions. Here we will 
focus on the trace dilation 

(1.1) K(u,ft) = ||K u (x)|| L oo (n) with K u (x) - l ' / '"' r)| 



(N 



(det du(x))™ 

Other dilation functionals used in the literature are the outer, inner and linear dilation (see 
[36J for more details) as well as mean dilations for mappings with finite distortion (see [1 Oj ) . 

There are a variety of extremal mapping problems in the theory of qc mappings, in fact qc 
mappings were introduced in just such a context in [22] ■ Extremal problems usually involve 
two domains Q,Q' C MJ 1 , (or two Riemann surfaces) for which there exists a quasiconformal 
mapping / : Q — > Q', and ask for a quasiconformal map u : Q — > Q' which minimizes 
a dilation function in a given class of competitors. Such competitors are usually other 
quasiconformal mappings with same boundary data as / on a portion (or all) of dQ or in 
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the same homotopy class as the given map /. Existence and uniqueness of extremals depend 
strongly on the dilation function used. Typically, existence follows from compactness and 
lower-semicontinuity arguments applied to a particular dilation function, and uniqueness 
does not hold unless the class of competitors is suitably restricted (for instance to Teichmuller 
mappings^) . 

Quasiconformal extremal problem arose first in the work of Grotzsch in the late 1920 's 
and were later studied in the two dimensional case both for open sets and for Riemann 
surfaces, see for instance [35], [2], [23], [33] and references therein. A celebrated result of 
Teichmuller, which was subsequently proved using two very different methods by Ahlfors 
[2] and by Bers [12] , states that given any orientation preserving homeomorphism / : S — > 
S' between two closed Riemann surfaces of genus g > 1 there exists among all mapping 
homotopic to /, a unique extremal which minimizes the L°° norm of the complex dilationj 
K(f,S) = \\Kf\\Loo/ff\. Moreover, the extremal mapping is a Teichmuller map, real analytic 
except at isolated points and with constant dilation Kf = const. In [23], Hamilton studied 
the extremal problem with a boundary data constraint, and one of his results is a maximum 
principle of sorts stating that if / is extremal, then the maximum of its Beltrami coefficient 
in S is the same as the maximum on OS. 

In higher dimensions, the problem becomes even more difficult and the references in the 
literature more sparse. The extremality problem without imposing boundary conditions is 
studied in the landmark paper [20]. Existence and uniqueness for the analogue of Grotzsch 
problem in higher dimensions is established in [18] and a maximum principle for C 2 extremal 
qc mappings is proved in [S] . More recently, in [TU] , [§] and [TJ the study of extremal problems 
for mappings of finite distortion is carried out for LP norms (and more general means) of the 
dilation functions with p finite, rather than with the L°° norm. In the same vein, the paper 
[TTj examines extremal problems in the mean for dilation functions based on the modulus of 
families of curves. 

In the literature discussed above, the study of extremal problems for qc mappings in space 
rests on a careful analysis of compactness properties for families of qc mappings with a 
uniform bound on dilation and on techniques from geometric function theory to establish 
uniqueness. The finite distortion problem relies on techniques from direct methods of calculus 
of variations, in which the study of the functional itself, rather than its Euler-Lagrange 
equations, is used. This approach is only natural as the extremal problem is posed in the 
class of qc mappings, and so there should be no additional hypothesis concerning second 
order derivatives. With this approach, however, there is so little regularity that finding 
information about the structure of extremal mappings (let alone the uniqueness) has proven 
intractable thus far. In particular, there is a huge gap between the findings in the two 
dimensional setting vs. the higher dimensional theory. 

In the present work we propose an approach to the extremal problem that is motivated by 
two classic papers: One by Ahlfors [2] in which an LP approximation of the L°° distortion 
is used to solve the extremal problem in the setting of Riemann surfaces. The other is by 



Roughly speaking, a planar qc mapping / is Teichmuller if there exist local conformal transformations 
(j>, ip such that 4> o / o ip^ 1 is afhne and <f> an d ip give rise to well defined quadratic differentials. 
2 The dilation Kt = |M ± M1 . 
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Aronsson [6], (see also [7]) where he assumes the extra hypothesis of C 2 regularity and carries 
out his program to determine the structure of absolute minimizing Lipschitz extensions. 

The extremal problem for qc mappings is a L°° variational problem that can be rephrased 
as follows: Given the boundary restriction Uq : dQ — > M. n of a C 1 (f2, MJ 1 ) qc mapping, find the 
qc extensions of Uq to Q with minimal trace dilation. From this viewpoint the problem has 
a superficial similarity to the problem of finding and studying minimal Lipschitz extensions 
u G Lip(Q) for scalar valued functions u G Lip(T) to a neighborhood r C O in such a 
way that Lip(u, Q) = Lip(uo, r)H The existence of minimal Lipschitz extensions was settled 
in 1934 by McShane (see also [15] for a more recent outlook of the problem), but simple 
examples show that uniqueness fails. In 1967, Aronsson showed that if the extremal condition 
is suitably localized to absolute minimal Lipschitz extension (AMLE), i.e., u G Lip(Q) is 
AMLE with respect to u G Lip(dQ) if Lip(u, V) = Lip(u, dV) for all V C Q, then a C 2 
function u is AMLE if and only if it solves the oo— Laplacian 

(1.2) UiUjUij = in Q. 

In essence, this PDE tells us that |Vw| is constant along the flow lines of Vtt. Aronsson 
also discovered several links between the geometry of the flow lines and the regularity and 
rigidity properties for oo— harmonic functions in planar regions. In the 1960's, solutions of 
( II. 2p could only be meaningfully defined as C 2 smooth. In the 1980's, however, a number 
of authors (see for instance [H], [25]) developed the theory of viscosity solutions, leading to 
Jensen's uniqueness theorem for AMLE and for the Dirichlet problem for the oo— Laplacian. 
Recent, exciting extensions of Aronsson's work to the vector-valued case provide further 
links with qc extremal problems (see Sheffield and Smart's preprint [33J ) but, as the theory 
of viscosity solutions has no vector valued counterpart, the standing C 2 hypothesis is present 
even in these very recent developments. 

The similarities with the AMLE theory prompted us to study a local form of the classical 
extremality condition, in which the qc mapping is required to have minimum dilation in every 
subset of the domain with respect to competitors having the same boundary values on that 
subset. Our goal is to find an operator that plays an analogous role to that of the oo-Laplacian 
in the characterization of extremals and would provide a platform for the qualitative study 
of these mappings. The non-linear relation between the dilation of a diffeomorphism and the 
dilation of its trace on a hypersurface introduce further complications in our work. 

In order to be more specific about our results we need to introduce some basic definitions: 
If <f) is a n x n matrix of C 1 functions, then the Ahlfors operator S(<p) is given by 

(1.3) 5(0) = £±^-itr(0I 

(see [3], [29], [1]). If u : Q — > Q' is a C 1 (H) orientation-preserving diffeomorphism then it is 
quasiconformal and det du > e > 0. For such a mapping we define the normalized pull back 
of the Euclidean metric under u~ l as the Riemannian metric g~ x . In coordinates, the metric 



3 We have set Lip{u, Q) = sup x yen x ^ y K |2_^ fa)l . 
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is expressed by the matrix^ 



i.) 



The inverse 



(dudu T )ij du,ikdujk 

9ij ~ 



(det du) 2 l n (det du) 2 l n ' 

In [21] the metric g is called t/ie distortion tensor. As in the work of Ahlfors [2] we consider 
LP approximations 

inf [ E? p (x)dx, 

V 



n 



of the L°° variational problem (these approximations have been studied in depth in |10j). 
Let Q C W l be a bounded open set. An orientation preserving QC mapping u : Q — > 
MP is p— extremal if ||K u || iP (n) < ||lK?j||,Lp(n) for all orientation preserving QC mappings 
v : Q — > M. n with u = v on dQ. It is straightforward to derive Euler-Lagrange equations 
for the LP variational problem: Every orientation preserving p— extremal diffeomorphism 
u = (u 1 , ...,u n ) G C 2 (Q,M. n ) satisfies the fully nonlinear system of PDE 

(L p uY = npd 1 [K n u ^ 2 (S(g)du^ T ) ij ] =npd j \^ 2 S{g) ti dv? i \ =0 

in Q, for % = 1, ...,n. Here (du)^ denotes the ij entry of du~ l , g iJ is defined by (11.4j) and 
5(g) by ( 11. 3p . For C 2 smooth mappings with non-singular Jacobian, the operator L p can 
be expressed in the non- divergence form (L p u) 1 = A l h(du)Uj£. The quasi-convexity of the LP 
variational functional [21] implies that the system satisfies the Legendre-Hadamard ellipticity 
conditions (see Lemma 13.11) . Motivated by the work of Aronsson, we consider the formal 
limit as p — > oo of the PDE L p u = and obtain 

(1-5) = T f^l(S(g)du- 1 ' T ) ij d x .K u = 0, 

or equivalently S(g)VK u = 0, where g = ^^ /n (see Section H below). This PDE tells us 

that the trace dilation K u is constant along the flow lines of the rows of the matrix S(g)du~ 1,T 
(and their linear combinations with C 1 coefficients). Since the derivation of ( II .5p is formal, 
a priori there need not be any links between solutions of this PDE and the extremal problem 
for qc mappings. However, such links exist and are addressed by the main results of the 
present paper. 

Theorem 1.1. Let Q C M. n be an open set. Lf u G C 2 (Q,M. n ) is an orientation preserving 
diffeomorphism solution of L^u = in Q, then for any bounded subdomain D dVt, 

K(u,D) < supK M . 

dD 



This metric has the following property: for all V, W € T^^R" we observe that (V, W) g -i ( u ( x )) = 
^ dU (dItdu~ ll Y' > ^ ' H ence u '■ i^'dx 2 ) — !> (^(fi),^ 1 ) is a conformal map in the sense that (duV, duW) g -i = 
(detdu) 2 / n (V,W) E u C \. 



Moreover, ifn>3 and K n has a strict maximum on dD in the sense that K„(z) < sup 9£) K u 
for z G D, then 

— n — l n— 1 

(1.6) K(u, D) = sup K u < y/E(n - l)sr sup K u n dD , 

dD dD 

where ^ U ,8D denotes the dilation of the trace of u on dD (see Definition \6. 

Corollary 1.2. Given the hypothesis of the previous theorem, 

(1) if mm xem K u (x) > y/n, then 

minK u (x) = minK u (x); 

(2) 7/K u is constant with K M > yfn on dVt then K u is constant in Q. Moreover, if n = 2 
and u is affine and is not conformal on dfl, then u is an ajfine map. 

Theorem 1.3. If u G C 2 (Q,M. n ) nC 1 (f2,]R n ) is an orientation preserving diffeomorphism, 
such that for every D dVt and v G C 2 (D, IR n ) R C l (D, IR n ) orientation preserving diffeomor- 
phism with u = v on dD we have K(u, D) < ~K(v , 5) then L^u = in Q. If n > 3 and for 
every D C Q, 

— 1 71 — 1 

(1.7) K(u, D) < n"5= sup K u " dD , 

dD 

then L^u = in Q. 

Corollary 1.4. Letu,v G C 2 (D,W l )r\C l (D,W n ) be orientation preserving diffeomorphisms, 
such that u = v on dD. If L^u = L^v = in D then K.(u,D) = K(t>,Z)). 

These results echo some of the n = 2 theory, in particular the maximum principle for the 
dilation in Theorem 11.11 recalls Hamilton's result [23l Corollary 2]. The fact that the dilation 
is constant along flow lines of a conformally invariant set of vectors recalls the analogous 
planar result about dilation being constant along the image of lines under the action of the 
conformal mappings associated to the quadratic differentials of Teichmuller mappings (see 
[341 Page 175] for a more detailed description). 

Remark 1.5. Theorem 1 1 . 1 1 and (16. ip tell us that if L^u = in Q C M™, n > 3, then for every 
D C Q for which K u has a strict maximum on dD, u is a quasi-minimizer for the extremal 
problem for the trace dilation in D. In fact, if v G C 2 (-D,IR n ) R C 1 (Z),]R n ) orientation 
preserving diffeomorphism with u = v on dD, 

K(u,D) < ^{n-iy^ supK^, = ^(n-iy^ supK^, < y/n(n-l)~^n^K(v, D). 

dD dD 

On the other hand, Theorem 11.31 tells us that those diffeomorphisms that are minimizers 
for the extremal problem for the trace dilation on every subset D C Q are also solutions of 
LoqU = 0. This lack of symmetry in our result follows from the fact that the constants in 
(11. 6p and (II. 7p are different. While the constant in (11.61) seems to be sharp, we are confident 
that is possible to improve on the constant in (jl.7p and conjecture: If u G C 2 (Q,M. n ) then 
the condition L^u = in Q is equivalent to minimizing the dilation K(n, D) < ¥L(v,D), on 
any subset D CC Q, among competitors v G C 2 (D,M™) R C l (D, W 1 ) with v = u on dD. 
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A large class of solutions of L^u = is provided by observing that (in any dimension) 
the set of C 2 solutions of L^u = is invariant by transformations u = F o u and v = u o F 
with F conformal. In particular, all the known explicit extremal QC mappings (that we are 
aware of) have constant trace dilation and hence satisfy the PDE (II. 5p . 

Corollary 1.6. (1) Any Teichmiiller map of the form u := ip o v o (fi^ 1 with ip, conformal 
and v affine is a solution of L^u = 0. (2) the QC mappings u(x) = la;!"" 1 ^ for a > solve 
L^u = away from the origin. (3) Let < a < 2ir and (r,9,z) be cylindrical coordinates 
forx = (xi, ...,x n ) where x\ = rcos9, xi = rsin^ and Xj = Zj, 3 < j < n. The QC mapping 



solves L^u = away from the set r = 0. 

The proofs of Theorems 11.11 and 11.31 rest on the analysis of the flow lines of the rows of 
the distortion tensor S(g) and the geometric interpretation of L^u = 0. We show that if u 
is not conformal on the boundary then these flow lines fill (row by row) the open set. 

The smoothness assumptions we make here are not natural for the problem, as they do not 
guarantee the necessary compactness properties that we need to prove existence of extremals. 
However, in the spirit of Aronsson's work on C 2 AMLE, it is plausible that the study of C 2 
mappings can yield a measure of intuition for the general setting. 

We observe that in the proof of the first part of Theorem 11.11 the smoothness hypothesis 
can be decreased to W 2 ' p for p sufficiently high, using the work of DiPerna and Lions [TH] 
(see also [5]) on solutions of ODE with rough coefficients. In fact, we can rephrase the PDE 
( II. 5p in the following terms: A QC mapping u : Q — >■ Q' is a weak solution of L^u = in 
Q if the trace of the corresponding distortion tensor g is constant along flow lines of linear 
combinations of the rows of S(g). In this formulation, the components of du need only be 
in a suitable Sobolev space or in BV. At present we are unable to decrease the smoothness 
hypothesis to the natural category of QC mappings and still obtain the maximum principle. 

Although currently we do not know how to prove existence of solutions of L^u = or how 
to attack the extremal problems for a fixed homotopy class of qc mappings, we indicate a 
possible strategy for a proof which involves the construction of a competitor for u by flowing 
u along a gradient flow for the LP norm of the dilation, then letting p — > oo. The initial 
value problem we need to control is the following: 



where Q = Q x (0, T) and d par fl = Q x {0} U dQ x (0, T). We prove the following 

Proposition 1.7. Let uq : Q — > M. n be a C 2,a diffeomorphism, for some < a < 1 with 
det duo > e > inQ. Assume that 



(1.8) 




(1.9) 




A#(duo)a#ug = o, 



for all x G dQ and i = 1, . . . , n. 

For every ji G (0, a) there exists positive constants C > depending onp, n, Q, e, llwollc 1 '"^) 
andT > depending on p,n,Q,e, \\uo\\c 2 ' a (n) an d a diffeomorphism u G C 2 ' fl (Q) solving (II. 9p 
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such that 



(1-10) |H| C 2, M(Q) + \\d t u\\ c o,^Q) < C\\u \\ C 2, a{n) , 

(1.11) det du > | /or a/Z (x, t) e 

We remark that although flows of qc mappings have been studied and used several times 
in the literature, see for instance [3], [29], [I], [13], and [32], this is the first instance of 
a gradient flow used in this context. Study of this flow may also contribute to a better 
understanding of the well-posedness and long-time behavior of initial value problems related 
to gradient flows of quasi-convex (and non convex) functionals (see [T7]). 

Acknowledgments. It is a pleasure to thank Hans Martin Reimann and Jeffrey Rauch for 
their interest and encouragement for this project. L. C. would like to dedicate this paper in 
fond memory of Juha Heinonen, who continues to be an inspiring role model. 

2. Preliminaries 
A map F : M™ — > M™ is conformal if at every point 

dF T dF = A/ n , 

for some scalar function A. Liouville's theorem states that if n > 2 then 1— quasiconformal 
mappings are conformal and that the only conformal mappings are compositions of rotations, 
dilations, and the inversion x i— > x/\x\ 2 . If n — 2, then orientation preserving conformal 
mappings are biholomorphisms (and vice versa). A simple computation shows that A = 
\dF\ 2 /n and det dF = vA™- We now list some equivalent formulations of conformality. 

Lemma 2.1. Let F : MJ 1 — > W n be a diffeomorphism. The following are equivalent: 

(a) F is conformal; 

(b) K F = ^Jn identically; 

(c) The expression (gLF) ji — n ^ffjfi vanishes identically; 



\dF\* 

\ — n 

(det dF) 2 / n < 



(d) = 0. 



Note that if n = 2 and u is holomorphic with du/dz ^ 0, then (duY 1 — n ^| = is a 
restatement of the Cauchy-Riemann equations. 

The action of conformal mappings on S, K u and g follows immediately from the definitions. 

Lemma 2.2. Let u : M. n — >• W 1 be a diffeomorphism and F be an orientation preserving 
conformal mapping. If we set u = F o u and denote by K and g the corresponding dilation 
and distortion tensor, then 

(a) K = K u ; 

(b) ~g = \- l dFgdF T ; 

(c) S(g) = \- 1 dFS{g)dF T ; 

(d) (du-T - n] §p = -nK-\dF T )^S{g){du^) T . 

In a similar fashion we will be interested in compositions with conformal mappings from 
the right, i.e., u = u o F, for which we can show: 
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Lemma 2.3. Let u : W 1 — > R" be a diffeomorphism and F be an orientation preserving 
conformal mapping. If we set u = u o F and denote by K and g the corresponding dilation 
and distortion tensor, then 

(a) K = K U ; 

(b) g = g; 

(c) S(g) = S(g); 

(d) (du-y -nj§s = -nK- 2 S(g)(du-Y(dF T )-\ 



3. The Euler-Lagrange system 



Let n C R" be a bounded, smooth, open set and u : O — > M. n a smooth, orientation 
preserving diffeomorphism with < det du < oo. For 1 < p < oo, we define, whenever the 
expression if finite, 



I I *j f2 



For any ip G C^°(fl, R n ) we set h(s) := T v [u + sip, Q) and compute 
d 



ds 



h(s) 



s=0 



(3.1; 



i 



np(det du)- p \du\ np - 2 du ■ dtfj -pidetduy^d^icof du)ij\du\ np dx 



(det du)P^ [ h | det du\P 3 1 ^ 



where cof du denotes the cofactor matrix of du, so that (cof du) T du = det duL Define the 
operator L p on M n -valued functions by 



{L p u) 1 = - pdj 



(du) ji - n 



(du)jj 
\du\ 2 



\du\ 



up 



(det du)P 



—pdj I du 1 



n- 



dudu 
\du\- 



T 



\du\ 



up 



(det du) 1 



j i 



(3.2) 



npdj K™ p ~ 2 (du~ 1 S(g)) :j . = npd j K^' 2 S(g) u du^ 



where du 13 denotes the ij entry of the inverse of du, and /„ is the nx n identity matrix, and 
K u is defined in (II. ip . g %3 by (II. 4p and S(g) by (11.31) . Note that the equality of the first and 
third expressions in ( 13. 2\\ uses 



(3.3) 



(du' 1 ) 1 -n 



du 
\du\' 



-nKZ 2 S(g)(du 



-1\T 



We write (L p u) 1 = djAUdu) where 



A)(q) = -p 



q 3 -n-^r 
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up 



(det q)P ' 



defined for any non-singular nxn matrix q. Notice that A 1 - (q)qij = 0. Set A l j k e (q) := -T^Aj(q). 
Recalling that 

d qM (coi q)ij = cof q M q 31 - cof q ie q jk and d gke q 31 = -q h q jk , 
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we compute 
(3.4) 

I lnp-2 



np(q ki q ji + q l0 q lk ) - n{np - 2)^ - \q\ 2 {q H q jk + pq £k q ji ) - nS^Sj 



\Q 



For C 2 smooth maps with non-singular Jacobian, the operator L p can be expressed in non- 
divergence form: 

(3.5) (L p uY = A%{du)u%. 

We remark that, in this form, the operator satisfies a Legendre-Hadamard ellipticity con- 
dition. This result can be inferred by observing that the functional T p {u, fi) is quasi-convex 
(it is actually polyconvex, this is proved in [2U Corollary 8.8.1]), and consequently, given 
sufficient smoothness, satisfies Legendre-Hadamard conditions. As we need explicit expres- 
sions for the constants involved, we provide the following estimates, whose elementary proof 
we omit. 

Lemma 3.1. For n > 3 and p > 1 or n > 2 and p > 1 and for all non-singular matrices q 
and vectors £, n G M. n , we have 

\ n \np-2 / \ \np-2 I |n(p+2)-2 

(3.6) ^(n^pfoHef^^ 1/1 



(d^)?^W'M I'/l 1^1 ^(del^^ (detg)^ 

where we can choose Ci(n,p) = n for n > 4 and p > 1 and /or n > 3 and p > 1; Ci(n,p) = 
if n = 3 and p > 1 and C\(n,p) = 2^j- /or n = 2 and p > 1. Tne constant C^n) does 
noi depend on p and can be chosen to be C^n) = 100n 3 . 

Remark 3.2. The operator L p does not satisfy the stronger ellipticity condition A|n| 2 > 

AftVijVht >Mv\ 2 - 

As the dilation functional is invariant under the action of conformal mappings (i.e., 
J-p(u,f2) = J-p(.F(n), f2) for all conformal mappings F : R n — > M. n that map f2 into itself), 
we can expect a corresponding invariance for the solutions of L p u = 0. 

Proposition 3.3. Let u : Q — > M™ 6e a orientation preserving diffeomorphism. 
(i) // F : M n — > W 1 is a conformal map and u = F(u), then 



(L p uY = ([dF-XY L p u^j 



where 



(L p v,y = —pdj ( du 1 



I„ — n 



dudu 



r 



\du\ 



rip 



j 1 



\du\ 2 J (±detdn)P 

with the sign in the denominator being +1 if F is orientation preserving and —1 
otherwise. 

(ii) // F : Q — )■ a composition of dilations, translations, and the inversion x i— )■ x/|x| 2 7 
t/ien v = u o F satisfies 

(L pV y = (L p uf 
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Remark 3.4. Case (i) holds for all conformal mappings, including in n = 2 all invertible 
holomorphic and anti-holomorphic functions. In contrast, case (ii) only applies to the given 
set of conformal transformations, as in the plane it fails to hold except for linear invertible 
holomorphic and anti-holomorphic functions. 

4. The Aronsson-Euler-Lagrange system and the operator I ra 
In this section we assume that for each p > 1 we have a solution u p of the PDE 

(4.1) LpUp = 0, in fl, 

and that u p — > n^ in C 2 norm on subcompacts of Q. Our goal is to formally derive a system 
of PDE for Uoo . 

Observe that dj\du\ np = np\du\ np ~ 2 u\u k ^ and dj \ det du\~ p ~ l = — (p+l)| det dn| _p_2 (cof du)kti 
Using the fact that d 3 (cof du)ij — for i — 1, . . . , n, we compute 

(4.2) (LpuY = -p ^ tdu ^ l^ jet^L ( M ^ C ° f du ^ ij + U ^° oi du ^ H 

( \du\ 2 \ 2 1 
— (p + 1) ( - — — j (cof du)ki(coi du)ij — n(np — 2)u k u l j — n\du\ 2 8ij8ik >u k j 

\ Lie L CLLL / J 

for all % = 1, . . . , n. 

Dividing the expression above by p 2 and letting p — > oo, we obtain that equation 

(14. ip formally converges to 

where 



n 



du 



2 \ 2 



J — — ( (cof du)ijU k + (cof du)kiu % - J — ( -j — ^— j (cof du)ij(coi du)ki — n 2 u l jU k 

Civ; U dLL \ J \ Clc L Li LA 



(4.3) 

= (ndu-ij — {du^du^foduki — \du\ 2 du ek )djduki- 
Observe that the system does not satisfy the Legendre-Hadamard conditions. 
Proposition 4.1. (1) Let u(x) = \x\ a ~ 1 x where a 6 R and a^O. T/ien 

Loou(x) = 

and 

n + a 2 — 1\ ~n(a 2 — l)(n — 1) x 



Zvufx ) 



p v y V a 2 J (n + a 2 -l)a \x\<*+^ 
away from the origin. 

(2) If u(r, 9, z) is defined by U.8\) from Corollary 1 1.6] then L^u = in the set r / 0. 
Proof. For (1), direct computation yields K u (x) 2 = """j^" 1 for all x ^ and 

^,2 





_ V\ 


V M 2 


n y 



S (9)v = 

The proof follows from these identities and from the definition of I ro and L p . 

For (2), for the case < 9 < a, we have detdw = it /a and |dn| 2 = (n — 1) + 7i 2 /a 2 . The 
computation in the a < 9 < 2n case is similar. □ 
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Note that u(x) = \x\ a x is conformal exactly when a = ±1, the only cases for which 
L p u = 0. 

5. Extremal mappings and the equation L^u = 
In this section we establish some analogues of Aronsson's results in [SI Section 3]. 
Lemma 5.1. If u € C 2 (Q,R n ), then 

(5-1) (L^uY = r ^^(S(g)du- 1 ' T ) i .d Xj K u . 

Proof. Observe that 

(5.2) fi^K, = ~(n^% - du^K u = K~ l (S(g)du~ l,T ). j , 



n \ \du\' 



and d Xj 'K u = (d qM K u )Uj e . The result follows quickly from ( 13 . 3H and ( 14. 3ft . □ 

The following proposition on conformal invariance of solutions of L^u follows immediately 
from combining previous lemma with Lemma 12.21 and Lemma 12.31 

Proposition 5.2. The set of C 2 solutions of L^u = is invariant by transformations 
u = F o u and v = u o F with F conformal. 

Corollary 5.3. In the plane any Teichmiiller map of the form u := ip o v o with ip,<j) 
conformal and v affine is a solution of L^u = 0. 

Lemma 5.4. IfK u = Kq > yfn, then there exists e = e(fTo) > so that 

e<\S(g)\ 2 <K(l-±-). 

Proof. Let < Ai < • • • < A n be the eigenvalues of g. We can write S(g) — g — —I and 
K 2 = tr(g). Direct computation yields 



trg 



\S(g)\ 2 = tr( g--2-I = tr(g 2 ) - - tr(g) 2 . 



71 



2\ 1 



71 



Note that the upper bound for \S(g)\ 2 is now immediate. 

We now prove the lower bound. Note that the n-tuple of positive numbers (Ai, . . . , A n ) 
satisfy Ai • • • A„ = 1 and K 2 = Ai H — ■ A n > n with equality if and only if Ai = ■ ■ • = A n = 1. 
Set A = K 2 Jn = ±(Ai + • • • A n ). Since tr{g 2 ) = J27=i \ 2 > it follows that 

£ (A, - A) 2 = £ \ 2 " 2A $> + n\ 2 = tr(g 2 ) - i tr(g) 2 = \S(g)\ 2 . 

i=l i=l i=l 

We now claim that for S > 0, there exists e = e(S) so that whenever K 2 /n = A > 1 + S then 
EA-A) 2 >e. 

To prove the claim, we argue by contradiction. Assume that there exists So > such that 
for each k G N we can find positive Af as in the hypothesis with X k — 1 > So > and 

(5.3) 

i 
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If A fc is a bounded sequence then so are Af (as Af > 0), hence for an appropriate subse- 
quence we may assume that A fc — > A > 1 + 5q and Af — > Aj as k — > oo. As Aj • • • A^ = 1 
for all it follows that Ai • • • A n = 1 and Aj > 0. From (15. 3 p we conclude that Aj = A and 
1 = Ai ■ ■ • A n = A n > (1 + So) n , a contradiction. 

If A fc is an unbounded sequence then for each M > there exists I — £m > such 
that A^ > M. On the other hand, in view of (I5.3P we have \\ > M/2 and consequently 
1 = \[ . . . \ e n > (M/2) n , a contradiction. □ 

Remark 5.5. When n = 2, we can find an explicit lower bound. In this case, A1A2 = 1 and 

\S(g)\ 2 = Xi + X 2 2 - ~(A X + A 2 ) 2 = i(A x + A 2 ) 2 - 2AiA 2 = ^(K* - 4). 

We are now ready to study the relation between C 2 extremal quasiconformal mappings 
and the operator L^. 

Proposition 5.6. Let Q C W 1 be an open set. If u G C 2 (f2,lR n ) is an orientation preserving 
diffeomorphism solution of L^u = in Q then for any bounded sub-domain D C Q, 

sup K u < sup K u . 

D 3D 

Proof. Let = sup 9D K u and assume that there exists po G D such that K u (po) = ko > ^ > 
y/n. Since u G C 2 and det du > e > 0, S(g)(du~ l ) T is Lipschitz in D. Consequently, for each 
Po £ D and i = 1, there exists a unique trajectory 7i(s) defined for s G / C 1 through 
Po satisfying (s) = [5(g) (dM^ 1 ) 7 ]^-^^)) f° r J — 1, Using (15. ip and the fact that 
LooM = 0, we have 

£*.(•*(.» = £ 

so 

K tt (7 i (s)) = Ku(po) 

for all s G I and all z = 1, . . . , n. If a curve 7$ terminates at a point p inside D, than at 
p there must exist another flow curve 7; that flows out of it. In fact, not all 7$ can have 
vanishing speed simultaneously at a point inside D. Arguing by contradiction, if this were 
to happen then we would have S(g) = at the end point. This would would yield K u = ^fn 
at the end point, while K M (7j(s)) = k Q > y/n, a contradiction. We choose i so that 



-Kn(7*(*)) = -J(^ Xj ^u(l(s)) = Sig^du-^KMs)) = 0, 



sup 

3 



S(gh(%(s)) >C n \S(g)\>0, 



for C n > 

The argument yields a piecewise C 1 curve 7 inside D, passing through p with K u (j(s)) 
k with 

^-iKs) = [S{g){du~Y]Ms)) 
for some index % = 1 , . . . , n and 



(5.4) sup S(g)Ms)) >C n \S(g) 



3 



for all s G /. There are two alternatives: (i) the curve 7 has finite length and so touches 
the boundary 3D in two points P, Q G D; (ii) The curve 7 does not touch dD and so has 
infinite length. 
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In (i), it follows that K U (P) = k > sup^K^ > K U (P), a contradiction that ko > /i. 

We need to exclude the second alternative. For simplicity we assume that the composition 
of flow lines is actually one single flow line, the general case is proved in the same way. For 
each i = 1 n, we have 



ti*( 7 (t)) - 
(for some I = 



d_ 

ds 



u\^(s))ds 



,71 



S(g)(du 



-1\T 



(j(s))duij(j(s))ds 



Hi 



(5.5) 

Consequently, for some > t[ > t, 



S(g)nh(s))ds 



sup 
=1,...,) 



K T (t))- u > )| >t 



sup 

= 1 ! 



\S(g)u\(7(ti)), 



and by (15. 4p . we conclude 



sup 

=!,...,) 



|u 4 (7(*))-« 4 (p«>)| >C n \S(g)Mt t ))t. 



Since ^(sOI is bounded from below by Lemma EH l M (7(^)) — M (Po)| has at least linear growth. 
Consequently, if 7 has infinite length, then u(D) would have to be unbounded, whereas since 
D is bounded so is u(D). □ 

We can now prove Corollary 11.21 

Proof of Corollary II. £1 Using the argument from the previous proof, we have that the set of 
points x G Q with K u (s) > ^Jn can be covered by compositions of flow lines of the rows of 
S(g)(du~ l ) T with K u constant along these curves. We have shown that if K u > ^Jn on such 
a curve then it must reach the boundary dVt. To prove (1) we observe that for any e > 
such that K u > y/n + e on dfl, if xq G {x G fi|K u (x) G (y/n, y/n + e)} then there exists 
a composition of flow lines passing through Xq which must reach the boundary and hence 
contradict the hypothesis K u > ^/n + e on dVL. As for (2), we observe that by virtue of (1) 
every point in Q can be connected to the boundary with a composition of flow lines along 
which K u is constant, thus concluding the proof. □ 

Remark 5.7. Arguing as in the proof of (15. ip . we can show that for each i = 1, .. 
let 7 : [0, e) — > Q be a flow line of the i-th row of S(g)du~ l,T , then for any j = 1, 
< t < e we have 

duij^it)) - ^(7(0)) = / —u)(-f(s))ds = j ^ k u) k (-i(s))ds 



n, if we 
.., n and 



(S(g)du- 1 ' T ) tk u} k ( 1 (s))ds 



K u d x K u (^(s))ds. 



This formula allows us to recover the differential of u from the dilation and the flow lines of 
the distortion tensor. In particular, if K u is constant in Q then the rows of du are constant 
along the flow lines of the corresponding rows of S{g)du~ 1 ' T . 



The previous remark yields: 
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Proposition 5.8. In the hypothesis of the previous theorem, if Q C M 2 and du (and hence 
K uy ) is constant in dfl, with K u > yfn on dfl, then du is constant in Q and hence u is affine. 

Proof. The remark above implies that if WL U is constant then the rows of du are constant 
along the flow lines of the corresponding rows of S(g)du~ 1,T . It suffices then to show that 
for every point po G fl we can find flow lines of both rows of S{g)du~ 1 ' T passing through 
that point and touching the boundary dfl. To establish this fact we recall that \S(g)\ > 
in fl and that, since we are in the planar case, both rows of S(g) cannot vanish unless they 
vanish simultaneously, which is impossible. Since du is invertible the rows of S(g)du~ x,T 
cannot vanish at any point in Q. Repeating the argument in the proof of Theorem 11.11 we 
see that the flow lines of the two rows of S(g)du~ 1,T through po cannot end in fl, nor can 
they continue for an infinite time, hence they must reach the boundary in a finite time. □ 

Remark 5.9. If n = 2 and K u > y/n on dfl then L^u = actually implies that K u is constant 
along any path in Q. Hence, in the plane there will be no C 2 (Q,M. 2 ) fl C 1 (fi,IR 2 ) solutions 
of LooU = in Q unless Kulan = const. 

We conclude this section with the proof of the necessity of the condition L^u = for a 
C 2 qc mapping to locally minimize dilation in subsets D C fl, among competitors with the 
same dilation on dD. 

Proposition 5.10. Let u G C 2 (Q,M. n ) be an orientation preserving diffeomorphism which 
does not solve L^u = in a closed ball D C fl. There exists v G C 2 (D,M. n ) orientation 
preserving diffeomorphism with u = v on dD such that K(t>, D) < D). 

Proof. Let u G C 2 (Q,W n ) be an orientation preserving diffeomorphism which does not solve 
LooU = in a closed ball D C fl. In view of the conformal invariance of the PDE we can 
assume without loss of generality that D = 5(0,1). Let E = {x G D\K(u,D) = K u (x)}. 
Since V X K U = at any interior point x G E, we must have E C dD and consequently 
~K(u,D) = sup{K u (x)\x G dQ}. If n denotes the outer unit normal at x to dD, then the 
latter yields that V x K u (x) = an for some a > at each x G E. The identity (15. ip then 
implies 

(5.6) S(g)du- 1 ' T n ^ 0. 

For A G R and x C 2 (D, IR n ), vanishing on dD, we define U\(x) = u(x) + A%(x). Using 
(15. 2p and a Taylor expansion of K UA in A, we have that 

(5.7) K ux =K U + Xd^Kudxij + 0(X 2 ) = K u + AK; 1 (S(g)du- l ' T ). j d Xij + 0(A 2 ). 

We claim that given u satisfying (15. 6p . we can find a mapping \ G C 2 (D,~R n ), vanishing on 
dD, such that the coefficient of A in (15 ,7p is strictly negative in a neighborhood U of E, for 
small values of A. This fact would allow us to the conclude the proof of the proposition. 
Indeed, for x G U H D and small values of A, we would have K UA < K u < K(w, D). On the 
other hand, for x G D\U, there would exist e > such that K u < K(u, D) — e, thus yielding 
that K Ux < K(u, D)-e + C\< K{u, D) for small values of A and C = C(||m|| c i, \\x\\c^D)- 
Given such inequalities we would then conclude that v = u\ is a qc diffeomorphism with the 
same boundary data as u and strictly smaller dilation K(«a, D) < K(w, D). 

To find X-, observe that if p G E then as a consequence of (15. 6p there exists v G M™ such 
that 

(5.8) (S{g)du- 1 ' T n,v) > 

14 



in a neighborhood B(p,r). Since we can cover E with a finite set of such neighborhoods, 
we obtain vectors V\,...,Vk G R™ for which (15. 8 j) holds in B(pk,r) and such that C 
Uf=i ^(Pij r )- For each Z = l,...,/c, let 0; : S" 1 ^ 1 — )■ R be a positive smooth function such 
that = outside B(p { ,r) n S n_1 . We set 



(5.9) X (x) = (1 - |x 



2\ 



i=l 



.r 



Clearly this mapping vanishes on dD and it can be easily modified near the origin to yield 
a smooth mapping in D. Observe that at every point in S n ~ l , 



k 

d X = -2(^1^1) 
^ 1=1 ' 



n. 



Substituting the latter in (15. 7B we obtain that for every point in dD, 
(5.10) 

k 

K ux = K u -2XK-\S(g)du~ 1 ' T ) ij n J v l ^ l +0(X 2 ) = K u -2XS^ 1 ^2(S{g)dvr 1 ' T n,^)<fn+0{\ 2 ). 

1=1 

In view of (15. 8 p and the choices of (pi and vi, it follows that for all x in a neighborhood 
E PI B(pi, r) C B(pi, r) R and A sufficiently small, the coefficient of A is strictly negative 
as 

k 

-2"KZ 1 ^{S{g)du- 1 ' T n,v l )(j> l < 0. 
i=\ 

Thus, the strict inequality K UA < K u holds, whereas elsewhere in dD \ uf =1 B(pi, r) we have 
equality. □ 

6. Dilation of traces of diffeomorphisms 

Throughout this section Q C R" is an open set, n > 3, u G C 2 (Q,M. n ) is an orientation 
preserving diffeomorphism, M C Q and M' = u(M) are closed, C 1 hypersurfaces endowed 
with metrics induced by the Euclidean metric. For x G M, we denote by ei,...,e n _i an 
orthonormal basis of T X M and by Wi, u> n _i an orthonormal basis of T u t x \M' . We let Wq 
be the unit normal field to M' such that (dun, w ) > 0. We denote by 



U = u\ M 

the trace of u on M. For each x G M consider the (n — 1) x (n — 1) matrix d M U(x) = (dy) 



with dij = (duei,Wj) 2 . 
Definition 6.1. The tangential dilation of U = u\m at a point x G M is given by 

|d M f/| 



1 

-1 



[dettZ M «]s= 

If f G C 1 (f2,R n ) is an orientation preserving diffeomorphism with u = v on M then 
Km,m = K«,Af on M. The following lemma is probably well known but we give a short proof 
as we did not find it in the literature. 
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Lemma 6.2. For every x G M , the dilation 

2 

ffii\ tl?"2 ^^ w ^x \dun\ 2 (dun, w )~^ 

(6.1) \ M <»"- lK « : — — ■ 

[det du\ 

Proof. We consider the two orthonormal frames of W 1 given by 

{n,ei,...,e n _i} and {«j 0i Wi,...,w„-i} 

and observe that in these frames du(x), x G M, can be represented as a block matrix 

(cfom, w ) 
{dun,Wi} d M U 



du 

Consequently, 



\du\ 2 = \d M U\ 2 + \dun\ 2 and det du = (dun,w ) det d M U. 

The estimate (16.11) then follows from the latter and from recalling (dun, wq) < y/n\du\. □ 

Lemma 16.21 and Proposition 15.101 immediately yield 

Proposition 6.3. If u G C 2 (Q,~R n ) is an orientation preserving diffeomorphism that does 
not solve L^u = in a ball D C Q then 

1 re. — 1 _ 

n'^ supK u ^ D < K(u,D). 

dD 

Theorem 11.31 now follows from Propositions 15.101 and 16.31 

We now turn to the final step in the proof of Theorem ll.il In order to estimate the dilation 
of the extension of u\m in terms of the tangential dilation we need more information about 
the extension. 

Lemma 6.4. Let u be a solution of L^u = in a neighborhood of M. If x G M satisfies 
VK u (x) 7^ and n || VK tt (i), then 

rn 1 2n 

(6.2) —Ku^jx) = K 2 , M (x). 

Proof. We observe that L^u = at x is equivalent to 

(6.3) du T dun \du\ 2 n = 0, 

n 

at x. In particular, n is an eigenvector of du T du(x) with eigenvalue \du\ 2 /n. Representing 
du 7 du in the orthonormal frame n, ei, e n _i, with eigenvectors of du T du(x), tangent to 
M corresponding to eigenvalues X 2 , % = 1, n — 1, we have the diagonal matrix 

^0 ... 



(6.4) du T du{x) 



T^rM = ( ^ ... 

... 



o ... 



We remark that du T du\ TxM = d M U T d M U, so that \d M U\ 2 = Ya=i x l From (E3D, we 
immediately obtain 

\du\ 2 = tr du T du = \^l + J2x 2 = l -^ + \d M U\ 2 



n * — ' n 

i=i 
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and 



det du 2 = det(du T du) = l^(det d M U) 2 . 



n 

To conclude, we have 

12 \du\ 2 



Km = — -^\du\—n-— = 1 - - )n-— 1 1 



(detrf-u)™- 1 V n / (detenu)" 



□ 



Remark 6.5. It is interesting to compare these conclusions with the example u(x) = \x\ a ~ 1 x 
on dB(0, 1). In this case, K 2 = n ~J 2 j'° and K„&b(oi) — n — 1. Note that the proof above 
does not apply as VK U = 0. 

Proof of Theorem The first part of Theorem 11.11 is proven in Proposition 15.61 For the 
second statement, observe that if x G 3D satisfies K(m, D) = sup dD K u = K u (x) then either 
VK u (x) = or it must be normal to dD. If V x K u (x) = 0, then the point x must be a local 
maximum of K u in Q. Consequently, there must exist a continuum F through x on which 
K u is constant and with FflD^O, otherwise x would be an isolated strict maximum point, 
an impossibility by the first part of Theorem 11.11 However the existence of points in D for 
which K u = sup dD K u contradicts the hypothesis K u (z) < sup^K^ for z G D and hence 
VK„(i) ^ 0. The proof now follows immediately from Proposition 15.61 and from (16.21) . □ 

7. QUASICONFORMAL GRADIENT FLOWS 

For a fixed diffeomorphism Uq : Q — > R ra , we want to study diffeomorphism solutions u(x,t) 
of the initial value problem (jl.9p . If there is a T > such that a solution u G C 2 (Q x (0, T)) 
exists with det du > in Q x (0,T), then by the same computations as in (13. ip . 

^F p (u,Q) = - ( j^r y |L p M| 2 cfe^ < 0, 

meaning that the p-distortion is nonincreasing along the flow. Hence we obtain 

Proposition 7.1. If u G C 2 (Q x [0, T), R n ) fl C X {Q, x [0,T),R n ) is a solution of flOD u»'to 
det du > zn fix [0, T), then for all < t < T, ||K Up \\ P LP ^ = \\~K. U \\ P LP ^— J Q \\L p u(-, t)\\L2(n)dt 
and consequently 

(7.1) ||l&u||LP(fix{i}) — W^uq \\lp(Q)- 

By Lemma 12.21 and Lemma 12.31 the functional T v {u, Q) is invariant by conformal defor- 
mation. Therefore, if we let s h-» F s : R™ — > R™ be a one-parameter semi-group of conformal 
transformations, then solutions to the PDE system 

d t u = L p u + ^-F s {u) 

would also satisfy (17. ip . Recall that the flow F s is conformal if 

S {dV) = dV + ^ - - trace (dV)I n 
2 n 
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where V = (£F S 



° ^T 1 = {£?.) 



s=0 



and S denotes the Ahlfors operator. If n = 2 then 



s=0 



this amounts to d{D = 0. If n > 3 there is more rigidity and conformality requires that 

T>(x) = a + Bx + 2(c • x)x — \x\ 2 c 



for any vectors a, c and matrix B with S'(-B) = (see 

We observe that in light of Proposition I3.3[ if u(x, t) is a solution of (jl.9p in f2 x (0, T) 
and v(x,t) = 5u(Xx,5~ 2 t) for some A, 5 > 0, then v(x,t) is still a solution with initial data 
Vo(^) = Suo(Xx) in an appropriately scaled domain. Applying inversions in a suitable way 
will also yield new solutions from u(x,t). 

7.1. Short-time existence from smooth initial diffeomorphisms. Throughout this 
section Q C M n is a bounded, C 2,a smooth, open set. 

Definition 7.2. Let H C I" be a smooth bounded domain and for T > let Q = Q x (0, T). 

The parabolic boundary is defined by d par Q = (Q x {t = 0}) U (<9f2 x (0,T)). The parabolic 
distance is d((x,t), (y, s)) := max(|x — y\, \J\t — s\). For a G (0, 1) we define the parabolic 
Hdlder class C°' a {Q) := {u G C(Q,R)\ \\ u \\c a (Q) '■— [ u }a + || M l|o < °°}) where 

\u(x,t) - u(y,s)\ 

(x,t),(y,s)eQ and (x,t)^(y,s) d((x, t) , (y, s)) a 

and |u| = sup Q \u\. For AT G N we let C K ' a (Q) = {u : Q ->• R| d Xi ■ ■ ■ d Xi u G C°' a (Q)}. 



[U\a '■ = SUp 



Proposition 7.3. Lei u : — ?■ M n fre a C 2,a diffeomorphism for some < a < 1 wrat/i 
det ciwo > e > inO. Assume that L p u = /or all x G <9f2. There exist constants C, T > 

Q = ft x (0,T) so tfia* 



depending on p,n,Q,e, ||«o||(7i><*ffi)> anc ^ a sequence of diffeomorphisms {u h } in C 2 ' a (Q) with 



(a) detu h > | for all (x,t) G 

(b) \\u h \\ C 2, a{Q) + \\d t u h \\ c o, a{Q) < C\\u \\ C 2, a{n) , 

d t u h * - Af l (du h - 1 )d j diu h ' k = in Q 
u h = Uq on d par Q. 



Proof. We prove the result by induction. We start with the base case h — 0. Since uo G 
C 2 ' a (Q), if we set a°f fc (a;) := A^f(c?u (a;)) then a°f fc G C 1,a (f2) and in view of Lemma |3TT| 

a°-'i k satisfies for all £, 77 G IR n and iffi 

(7.2) A*iei 2 M 2 < ojf tt ^i» < A h ierM 2 , 

with h = and 

, 7 ,x A A r , x 2( \du h \ np - 2 , IjUftT^! ^^/ x 2 (||^o||^(n))" (P+2) - 2 

(7 ' 3) = ^ l(det^ + (det^)^ J " CMP ^ 

and 

^h ^ , ^ \du h \ n P~ 2 ^ C{n) 



(detdM /l )P C/i||dwo|U«>(r2) 



with h = and Co = 1. We have also used the bound det q < n\q\ 
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Applying Lemma IA.1I with T = 1 we obtain a constant Co = Co(n,p,e, \\duo\\c a (n)) > 
and a map u 1 G C 2,a {Q) that solves (c) and satisfies 

(7.5) T a / 2 (T- l / 2 \\du h \\ ca{Q) + \\u h \\ C 2, a(Q) + || fyi fc ||co,« w) ) < C A _i||uo||cv. ( n), 
for /i = 1. 

Next, the bound on T will be imposed to keep the determinant from vanishing. Set 



w = u 



(7.6) 



Uo and observe that this map solves the equation 

d t w i - Af l (du°)d j diw i = Af^du^djdiui in Q 

w = on dparQ 



for i = 1, . . . , n. An application of the Schauder estimates flA.4j) yields 

1 — a 

(7.7) T —\\dw\\ c «{Q) < C(n,p,e, ||dito||c«(n))||%||c*.«*(n)- 

Choose T\ < 1 sufficiently small depending only on n, p, a, \\du \\c^(n) and e = min^ det du 
so that ||du>||c7«(Q) <C §■ Since the determinant has polynomial dependence on the coeffi- 
cients, we have (a) for h — 1 in Qi — $7 x (0,Ti). 

Next we iterate this process to generate u h from u" ,_1 , h > 2, yielding estimates of the 
form (17.51) in Q h = Q x (0,T^) for some constants Ch,T h > 0. The difficulty resides in 
controlling the constants Ch and Th independently of h. In the following lemma, we show 
how to (re) choose the constants Ch = £ and Th = T uniformly in h e N and while keeping 
det du h+l > e/2 in Q h = Q x (0, f ). 

Lemma 7.4. Using the notation of Proposition \773\ if there exist C, T > with T < 1 sttc/i 



(7.8) Hdw' 1 - dwo||c«(a) < e anc? ll^^^llc^cs) < C||u ||c 2 .«(n) 

/or /i = 1, N — l and Q = Q x [0, T], then there exist constants <£ = £(n,p, \\u \\ci,<*m)) 
and T > X = %(C,n,p, ||mo||c* 1 ' q (q)) > that are independent of N and such that 

\\du N -du \\ o <*(Q) < e and \\djdiu N \\c"(Q) < ^||mo||c 2 ^(Q) 

in Q = O x [0,2]. 



> 



Proof. We set 
(7.9) 



u iv _ M ^ i an( ^ Q^ggj-yg thai u) N satisfies 



d t w N ' { - Af^du^djdiw^ 



w 



N 







ik/J n ,N-2\ 



Af^du*- 1 ) - Aftidu 



d j diu N ~ 1 ' k in Q 



on <9 par Q. 



Applying the Schauder estimates flA.41) in the cylinder Q x [0, T] with < X < T to be 
chosen, we obtain 



\dw N \\c~ ( Q) < C(n,P,e, \\du N - l \\ Ca{Q) )%^ 2 



A%{du N - l )-A%{du N - 2 ) 



d j d l u N ~ 1 ' k 



C a (Q) 



The hypothesis (17.81) yields a bound on the Holder norm of the second derivatives 



C a {Q) < C||M ||c 2 .«(n)) 
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and at the same time a strictly positive lower bound on det du h > |, for h — 1, . . . , N — 1 in 
Q so that 

(7.10) H^llc^Q) < c(n, P ,e, ii^-Ic^q))^!!^!!^^)^ 1 -^ 2 !!^- 1 !!^^) 

< C(n,p,e, \\du \\c a (Q))C\\uo\\ C 2, a (Q)1 {1 ' a)/ ' 2 \\dw N ' 1 \\c^(Q)- 

Choosing % sufficiently small depending only on C,n,p and ||'Uo||c ,2 > Q (n) an d independent of 
the index N, it follows that 

IM^IIccq) < 9\\dw N ~ 1 \\c<*(Q), 
where 6 e (0, 1). The latter yields 

N 9 
\\du N - du \\c«(Q) < ^2 W du3 ~ du3 ~ l II C a (Q) - Y^6^ dul ~ du ^ Ca tQ)' 
i=i 

We have proved the first part of the conclusion. To establish the estimate \\djdiu N \\c<*(Q) < 
£||iio||c 2 >«(fi) it is now sufficient to apply Schauder estimates to (c) with h = N — 1 and 
observe that the ellipticity bounds on A and A are independent of TV in light of the estimate 
\\du N - du \\ c ^(Q) < e. □ 

We now complete the proof of Proposition [7[3j Applying Lemma [TT^l to the case h = N = 2 
yields bounds 

\\du 2 - du \\ C i(Q) < e and \\djdiu h \\c<*(cj) < C||uo||o 2 . a (n) 

in Q = Q x [0, T], with % = ^(Cilluollc 2 '"^); n -,V-, e) > and €, = C(n,p, e, H^oHc 1 -"^)) > 0. 
As £ is a constant independent of Ci, we can eliminate the dependence on h by applying 
Lemma [7.41 again, yielding 

(7.11) \\du h - du \\c<*(Q) < e and \\djdiu 2 \\c a (Q) < £|KHc 2 ."(c) 
for h = 2 in Q = Q x [0, f ] with 

f = f(€,n,p,e, ||wo||c 2 .«(c)) = T(n,p,e, \\u \\c^ {n) ) > 0. 

At this point we proceed by induction on h: If (17. lip holds for h — 1, . . . , N in Q — Q x [0, T] 
with T = T(n,p, e, || w 1| c ,2 - a (^) ) > and <£ = C(n,p, e, ||iio||cWn)) > 0, then applying Lemma 
17^1 at the level of + 1 leads to (I77TTD for h = N + 1 in Q = Q x [0, f ] with f and £ as 
above; there is no degeneracy of the constants. Finally, since % is uniform in h, (b) follows 
from the Schauder estimate fl A. T|) . This concludes the proof of the proposition. □ 

The previous proposition and Arzela-Ascoli theorem yields Proposition 11.71 

Remark 7.5. The proof of the short time existence is quite standard and uses only the 
Legendre-Hadamard ellipticity rather than the specific structure of the non-linearity in the 
PDE. It seems plausible to expect that techniques such as those in the paper may also 
be used in our setting to establish short-time existence for C l,a initial data. 

Note that the Schauder estimates in the appendix yield uniqueness of a C 2,fl solution (for 
short time) but nevertheless there still may exist rough solutions of the equations with the 
same initial data |28j . 
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Appendix A. Existence and basic estimates for classical solutions of 

PARABOLIC SYSTEMS 

We recall results of Schlag [31] and Misawa [27] concerning classical (i.e., two spatial and 
one time derivative in C a ) solutions of the systenu 



(A.l) 



d t w - Ajjftx, t)w^ = f{x, t) in Q 
w = on d par Q 



assuming; that VL is a C 2 ' a domain, Q = fix (0,T), A, f e C a (Q) } the compatibility condition 
/ = on dfl x {t = 0} and an ellipticity assumption 

(A.2) < A$(x,t)^ <A|£| 2 , 

for some A, A > and for all (x, t) G Q and £ G R nxn . 

Schlag proves that there exists a constant C = C(n,\,A,\\A\\o )(X ,Q) > such that if 
w G C 2 ' Q (Q) and d t w G C°' a (Q) solves then 

(A.3) [ Wjl ) a + [w t ] a <C(\w\ + [f) a ). 

Misawa proved that such solutions exist and that the estimate can be slightly strengthened 

ll^j/llc-CQ) + ll w t||c«(Q) + II Vw||c«(Q) + ||w||c«(Q) < C||/||c«(Q), 

with a constant C > depending on n, A, A, || A||c«(q), fi and T. 

We address the dependence of the constants in the Schauder estimates from the parameter 
T. Since these estimates have a local character we expect the constant to blow up as T — > oo 
and to be bounded for T > fixed. 

Let T > T > and set ^fi := {x G M n |VTj! G ft}. Observe that if w solves (TP) then 

the function w(x,t) := w(yTx,Tt) solves 

d t w - Af^Vfx, Tt)w k 3l = f(x, t) := Tf(VTx } Tt) in ^fi x (0, 1) 
w = on dpar-j^Q X (0, 1) 

Note that d t w(x, t) = Td t w(VTx, Tt), djiw(x, t) = Tdjiw(y/Tx, Tt) and Vw(x, t) = y/TVw(y/Tx, 
The Holder norm of A(x,t) := A(y/Tx,Tt) is bounded by 



o/2n 



min{l,r a / 2 }||A|| c , (0) < ||i|| c « ( ^ nx(0)1)) < \\A\\c<* iQ )(l + T j. 

Since the ellipticity constants of the coefficients are not affected by the rescaling, the Schauder 
estimates for w read 

||^i/||c«(^nx(o,i))+lNtllc«(^nx(o,i))+ll^^llc«(^nx(o,i))+ll^llc«(^nx(o,i)^ — C|l/llc a (^nx(o,i))' 

with a constant C > depending on n, A, A, ||A|| c^q), f2. Rescaling back this estimate to 
the parabolic cylinder Q = x (0,T), we obtain 
(A.4) 



min{l,T°/ 2 } 



(ikiillo-W) + IK||c«(Q) + T 1/2 \\Vw\\c«(q) +T l \\w\\c<*(Q) \ < C||/||c<*(q) 



with C depending on the quantities above and on T Q . 



5 Both papers address more general systems. 
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Using the standard method based on applying Fourier transform to the integral 

A i Jl(u k (j))i(u i (j)) j dx 

n 

(see for instance [2T]) we note that the Schauder estimates continue to hold when weakening 
the ellipticity assumption from ( 1A.2I) to the Legendre-Hadamard ellipticity 



(A.5) A|e| 2 |^| 2 < Af^xM^W < A \Z\ 2 \V\ 2 

for some A, A > and for all (x, t) £ Q and £, rj £ IR n . Recasting these results for the system 
(A.6) 



d t u % — A l J k l (x,t)d j diu k = in Q 
u{x, 0) = uq{x) for all x £ d par Q 



we obtain the following: 



Lemma A.l. Assume that d£l is C 1 ' , T > and for < T < T , Q = SI x (0,T). // 
A £ C°' a (Q) and the compatibility condition 

A % Jl(x, 0)djdiUQ(x) = for all x £ dQ and i — l,...,n, 

holds, then given u £ C 2,a (Q) there exists a solution u £ C 2 ' a (Q) of (1A.6|) withu t £ C°' a ((5) 
and such that 

( A - 7 ) ll M llc 2 .-(Q) + ||^«l|c«(o) < Ci||mq||c«(q)- 

The positive constant C\ depends only T,n,Q, \, A and the C 2 ' a norm of the coefficients of 

A. The time-scaled version of OA .70 is 

(A.8) 

min{l,T a / 2 }/ 1/2 x \ .. 

^75— (JK'/llc^Q) + IKIIc<*(Q) +T 1 \\\/u\\c<*(Q)+T ||u||c««(0)J < <^2|KI|c° 

1 + Tg 

where Ci depends only T , n, Q, A, A and the C a norm of the coefficients of A. 



*(<2)> 



Appendix B. Evolution equations for the Jacobian and the distortion 

TENSOR. 

Let u £ C 1 ([0, T], C 3 (Q, R™)) be a classical solution of ( II. 9p and Q be the range of u (or 
equivalently, the range of u(-,t) for all t £ [0,T]). Denote by v{-,t) = w _1 (-,t) the inverse of 
the diffeomorphism u at time t and set (3(y) = det dv(y). For a fixed time i, set y = u(x,t) 
and dv(y,t) = du~ l {x,t). Let £ £ Cg°(f2 X [0,T],R). 

The argument in [T?j Theorem 2.1] yields 

= / / [^it* - d x .A)(du)] {d Vi i)\ u dx dt 
Jo Jn 

-T 



(B.l) = / / [-dMu{x,t),t)) + d t u\d yi O\ u ]-d X] A){du){d yi O\ u dxdt 

Jo Jn 

= I I — d t £\ det dv\ — d Xj A l j(du) det dv\ {d yi £) | J detdiicfedt 

Jo Jn\ " u 3 / 
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Next, we define 



up 



4M = 4(^) = -p(*- i -«^ p )^ 



for all non-singular matrices q and observe that 

d x .A){du{x,t)) = d x .A){dv{u{x,t),t)) = dv h '(u(x,t),t)[d yh A i j (dv(y,t))]\ u . 
Also, on Q, 

\du\ np 



dv h >d yh A){dv)P = -d yh 



-a, 



Uh 



p(cof civ),-/, ( aV l - n 
p(cof (if ) J 7 i I aV* — n- 



„ (det 



-a 



<//, 



p/3 [Sfu-n 



duhjdui 



\du\ : 



\du\ 



v (det du)P 

np 



v (det du)P 



since c^ fc (cof (if )jh = 0. The latter and ( IB.lj) yield 



= ^ jf ( - I u det di; | u - dv hj \ u [d yh A) (dv (y , t ) )] | u det dv | U (9 W k) det rfM ^ 




- 



'0 ./n 




Uh 



o Jn 



P[5hi-n 



duhjduij 
\du\ 2 



d y £dy dt 
\du\ np 



(det du)P 



£ dydt. 



We have then proved the following: 



Lemma B.l. Letu G C 1 ([0, T], C 3 (f2, R n )) be a classical solution of §nfl). If we setv{-,t) 
u~ l (-,t) and /3(y) = det dv (y), then (3 satisfies 



(B.2) 

in Q x (0, T), with 



= dy.dy h 



B ih (du) 







B ih (du) = p 



duhjduij\ \du\ 
Ohi — n — i , | , \ ; 



np 



\du\ 2 J (det du)r 



as well as Neumann type conditions 



dV^dyMdv)]? 



d v d Vh 



B ih (du) 











for all (y,t) effix (0,T). 

Let rj e M. n and q a non-singular matrix, and consider the quantity 



|g| 2 /(detg) p 
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2 A _ A M 



ftp 



/ (det g)P 



In the model case when n = 2, q is diagonal with eigenvalues < Ai < A2, and for unit 77, 
one has 



B ih (q) V Y=P[l- 



Af+Aj 
2 



A1A2 



The matrix does not have a sign, it vanishes when Ai = A2. Unlike the case studied in 
|17j . the parabolic maximum principle does not apply. 

Lemma B.2. If u is as in Lemma \B.1\ then the corresponding conformal metric evolves 
according to 



(B.3) d t g aP = npb^ n \d k d j 



S{g) la dufK np - 2 



?? 



d Xj d Xk I S ih (g) detdu- l K np - 2 )du kh du jl -d Xk I S ih {g) det du' 1 K np ~ 2 )du sh du kl d Xj d Xs u l du 3 



dunk + d k dj 



S{g) ip dufK np -' 2 



du 



ak 



+ (det duY l du kl du is d k d s u l d j S(g) mi duj^K np - 2 



dUa^dUphx 



in Q — ft x (0, T) with g = g on d par Q, and where b = (det du) 1 = det du 1 o u. 
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